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Abstract. Given a nonnegative polynomial /, we provide an explicit expres- 
sion for its best £i-norm approximation by a sum of squares of given degree. 



1. Introduction 

This note is concerned with the cone of nonnegative polynomials and its sub- 
cone of polynomials that are sums of squares (s.o.s.). Understanding the difference 
between these two cones is of practical importance because if on the one hand 
nonnegative polynomials are ubiquitous, on the other hand s.o.s. polynomials are 
much easier to handle. For instance, and in contrast with nonnegative polynomials, 
checking whether a given polynomial is s.o.s. can be done efficiently by solving a 
semidefinite program, a powerful technique of convex optimization. 

A negative result by Blekherman [3] states that when the degree is fixed, there 
are much more nonnegative polynomials than sums of squares and the gap between 
the two corresponding cones increases with the number of variables. On the other 
hand, if the degree is allowed to vary, it has been known for some time that the 
cone of s.o.s. polynomials is dense (for the ^i-norm of coefficients) in the cone of 
polynomials nonnegative on the box [— 1, 1]". See e.g. Berg, Christensen and Ressel 
PQ and Berg [2]. However, [1] was essentially an existence result and subsequently, 
Lasserre and Netzer [5] have provided a very simple and explicit sequence of s.o.s. 
polynomials converging for the £i-norm to a given nonnegative polynomial /. 

In this note we provide an explicit expression for the best ^i-norm approximation 
of a given nonnegative polynomial / € R[x] by a s.o.s. polynomial g of given degree 
2d (> deg/). It turns out that 



9 = f + 



for some nonnegative vector A* € R n+1 , very much like the approximation already 
provided in [5] (where the A* 's are equal) . In addition, the vector A* is an optimal 
solution of an explicit semidefinite program, and so can be computed efficiently. 

2. Main result 

2.1. Notation and definitions. Let K[x] (resp. R[x]<j) denote the ring of real 
polynomials in the variables x = (xi, . . . ,x n ) (resp. polynomials of degree at 
most d), whereas £[x] (resp. E[x]<j) denotes its subset of sums of squares (s.o.s.) 
polynomials (resp. of s.o.s. of degree at most 2d). For every a 6 N™ the notation 
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x%> and for every i 6 N, let N p d := {/3 e N" : 
53j ft — ^} wno se cardinal is s(d) = A polynomial / € R[x] is written 

x->/(x) = £ / Q x Q , 

and / can be identified with its vector of coefficients f = (f a ) in the canonical basis 
(x Q ), a e N™. Hence, denote by ||/||i the £i-norm ^ Q \f a \ of the coefficient vector 
f which also defines a norm on R[x]d. 

Let S p C M. pxp denote the space of real p x p symmetric matrices. For any 
two matrices A,B S S p , the notation A b (resp. >- 0) stands for A is positive 
semidefinite (resp. positive definite), and the notation (A, B) stands for trace AB. 

Let Vd(x) = (x a ), a e N'J, and let B a € R s ( rf ) xs ( d ) be real symmetric matrices 
such that 

(2.1) v d (x)v rf (x) T = Y, X " B - 

Recall that a polynomial g € R[x]2d is a s.o.s. if and only if there exists a real 
positive semidefinite matrix X 6 R. S W XS W such that 

g a = (X,B Q ), VaeNfc. 

d-moment matrix. With a sequence y = (y a ), a £ N n , let L y : K[x] — > K be the 
linear functional 

/ (=^/ Q x Q ) ^ £y(/) = J^/alfa, /£%]■ 

With <i G N, the <i-moment matrix associated with y is the real symmetric matrix 
Mrf(y) with rows and columns indexed in N d , and defined by: 

(2.2) M d (y)(a,P) := L y (x a +^) - y a+p , Va.jSefl?. 
It is straightforward to check that 

{L y («? 2 )>0 V fl eR[x] d } M d (y) b 0, d = 0, 1, . . . . 

Semidefinite programming. A semidefinite program is a convex (more precisely 
convex conic) optimization problem of the form minx {(C, X) : „4X = b; X y 0}, 
for some real symmetric matrices C,X S S p , vector b g R m , and some linear 
mapping A : 5 P — » M m . Semidefinite programming is a powerful technique of 
convex optimization, ubiquitous in many areas. A semidefinite program can be 
solved efficiently and even in time polynomial in the input size of the problem, for 
fixed arbitrary precision. For more details the interested reader is referred to e.g. 

2.2. The result. Consider the following optimization problem: 

(2.3) pd-^miniWf-gWt : g e E[x] d }, 

9 

that is, one searches for the best t i-approximation of / by a s.o.s. polynomial of 
degree at most 2d (> deg /). Of course, and even though (|2.3j) is well defined for an 
arbitrary / G K[x], such a problem is of particular interest when / is nonnegative 
but not a s.o.s. 
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Theorem 2.1. Let f G R[x] and let 2d > degf. The best t\-norm approximation 
of f by a s.o.s. polynomial of degree at most 2d is given by 



(2.4) 



x -> ff (x) = /(x)+ (\* Q + J2\* 



„"2d\ 



for some nonnegative vector A* G R" +1 . Hence pd = A*, and in addition, X* 



i=0 



an optimal solution of the semidefinite program: 

( n n 

(2.5) min 



. i=0 



i=l 



Proof. Consider / as an element of K[x] 2 d by setting f a =0 whenever \a\ > degf 
(where \a\ = ^o^), and rewrite (|2 . 3[) as the semidefinite program: 



Pd 



mm 

A>0,XhO,c/ 



(2.6) 



s.t. X a + g a 
A Q — g a 



> f a , Va G 
>-/«, VaGN" d 



g a -(X,B a ) = 0, 
The dual semidefinite program of (|2.6p reads: 



Va G N 



2(i' 



max - v a ) 

u a ,v a >0,y — ' 
_ aGN™ 



S.t. 



u a + v a < 1 Va G N£ d 
M d (y) b 



VaeN^, 



(2.7) 

or, equivalently, 
(2.8) 



The semidefinite program (|2.8p has an optimal solution y* because the feasible set 
is compact. In addition, let y — (y a ) be the moment sequence of the measure 
dp = e~" x " dx, scaled so that \y a \ < 1 for all a G N^- Then (y, u,v) with 
u = — min[y,0] and v = max[y,0], is strictly feasible in (|2.7[) because M^y) >- 0, 
and so Slater's condition^ holds for (|2.7p . Therefore, by a standard duality result in 
convex optimization, there is no duality gap between (12 .6[) and (|2.7p (or (|2.8p ). and 
(|2.6p has an optimal solution (A*,X*,g*). Hence pd = —L y *(f) for any optimal 
solution y* of 



max —Ly(f) 

s.t. M d {y) t 

\y a \ < 1, VaGN" d . 



Slater's condition holds the conic optimization problem P : min x {c'x : Ax = b; x G K}, 
where K C 1" is a convex cone, if there exists a strictly feasible solution xo G intK. In this case, 
there is no duality gap between P and its dual P* : max z {b'z : c — A'z £ K*}. In addition, if 
the optimal value is bounded then P* has an optimal solution. 



4 



JEAN B. LASSERRE 



Now by [6l Lemma 1], M<j(y) h implies that \y a \ < max[L y (l), max^ L y {xj )], 



for every a G N 



2<l 



Therefore 

Pd : 



(2.9) 



mm 
y 
s.t. 



has the equivalent formulation 

Ly(f)) 

M d (y) >= 
L y (l) < 1 



1. 



and it has an 
e Rl +1 with 

□ 



whose dual is eaxctly (|2.5I) . Again Slater's condition holds for (|2.9[) 
optimal solution y*. Therefore (|2.5[) also has an optimal solution A" 
Pd = J2i K i tne desired result. 

So the best £i-norm s.o.s. approximation 5 in Theorem (|2.1I) is very much 
the same as the ^-approximation provided in Lasserre and Netzer [5] where all 
coefficients A* were identical. 

Example 1. Consider the Motzkin-like polynomial x n> /(x) = x\x\{x\ + x\ — 
f) + 1/27 of degree 6, which is nonnegative but not a s.o.s., and with a global 
minimum /* = attained at 4 global minimizers x* = (±(1/3) 1/2 , ±(1/3) 1/2 ). The 
results are displayed in Tabled] for d = 3,4,5. 



d 


A* 


Pd 


3 


ps 10- 


a (5.445, 5.367, 5.367) 


1.6 nr 2 


4 


ps 10' 


4 (2.4,9.36,9.36) 


PS 2.10" 3 


5 


w 10" 


5 (0.04, 4.34, 4.34) 


ps 8. 10- 5 



Table 1. Best ^-approximation for the Motzkin polynomial. 
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2 Computation was made by running the GloptiPoly software [4] dedicated to solving the gen- 
eralized problem of moments. 



